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Abstract 



lf~) \ We study rare exclusive B decays based on the quark-level transition b — ► 

[ s(d)l + l~, where I = e or /i, in the context of super symmetric theories with 

minimal flavour violation. We present analytic expressions for various mixing 
matrices in the presence of new CP- violating phases, and examine their impact 
on observables involving B and B decays. An estimate is obtained for CP- 
CN \ violating asymmetries in B — > K^*'l + l~ and B — * p(tt)1 + 1~ decays for the 

| dilepton invariant mass region 1.2 GeV < M\+i- < Mj/^. As a typical result, 

we find a CP-violating partial width asymmetry of about —6% (—5%) in the 
■ case of B — > ir (B — > p) in effective supersymmetry with phases of O(l), taking 

into account the measurement of the inclusive b — > S7 branching fraction. On 
Qh. the other hand, CP asymmetries of less than 1% are predicted in the case of 

rn ' B -» i^W. We argue that it is not sufficient to have additional CP phases of 

O(l) to observe large CP-violating effects in exclusive b — ► s(d)l + l~ decays. 
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I. INTRODUCTION 



Within the standard model (SM), CP violation is caused by a non-zero complex phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [|l|. While the experimentally 
observed indirect CP violation in the neutral kaon system, ex, can be accommodated in the 
SM, it is still an open question whether the SM description of CP violation is consistent with 
the new experimental result on direct CP violation, e'/ex, since the theoretical prediction 
of its precise value suffers from large hadronic uncertainties 0. On the other hand, if the 
baryon asymmetry of the universe has been generated via baryogenesis at the electroweak 
phase transition, the CKM mechanism of CP violation cannot account for the observed 
amount of baryon asymmetry. This feature could be a hint of the existence of CP- violating 
sources outside the CKM matrix |J. 

Important tests of the SM are provided by flavour-changing neutral current (FCNC) 
reactions involving B decays ||, thus offering an opportunity to search for supersymmetric 
extensions of the SM |5|,|6[]. There are at present only a few FCNC processes which have been 
observed experimentally, but the situation will change considerably after the completion of 
B factories in the near future. 

In this work, we analyse the exclusive decays B — > K^*' l + l~ and B — > p(ix) l + l~ in 
the context of supersymmetry (SUSY) with minimal particle content and R-parity conser- 
vation The inclusive reaction B — > X s l + l~ within supersymmetric models has been 
extensively studied in Refs. and, more recently, in Ref. pTTfl . New physics effects in 



the exclusive channels have been investigated in Refs. [jT8| pOf - 

We place particular emphasis on CP-violating effects associated with the partial rate 
asymmetry between B and B decays as well as the forward-backward asymmetry of the 
l~ . Within the SM these effects turn out to be unobservably small 10~ 3 ) in the decays 
B — > K^l + l~ plj, and amount to only a few per cent in B — > p(n)l + l~ |2T|, p2"f . However, 



in models with new CP- violating phases in addition to the single phase of the CKM matrix, 
larger effects may occur due to the interference of amplitudes with different phases. The 
purpose of the present analysis is to explore CP-violating observables in the aforementioned 
FCNC reactions that could provide evidence of a non-standard source of CP violation, and 
hence may be useful in analysing supersymmetry in future collider experiments. 

The paper is organized as follows. In Sec. O, we exhibit the various mixing matrices of the 
minimal supersymmetric standard model (MSSM) in the presence of additional CP-violating 
phases. Within such a framework we discuss different scenarios for the SUSY parameters. 
In Sec. [HI] , we are primarily concerned with the short- distance matrix element and Wilson 
coefficients governing b —>■ s(d) in the MSSM. We also briefly describe an approximate 
procedure to incorporate quark antiquark resonant intermediate states - namely p, u, and 
the J/ip family - which enter through the decay chain b — > s(d)V Q q — > s(d)l + l~. Section 
[P7| is devoted to the exclusive decay modes B — > K(*H + l~ and B — > n(p)l + l~, where 
formulae are given to calculate CP asymmetries which can be determined experimentally by 
measuring the difference of B and B events. In Sec. [V], we present our numerical results for 
CP-violating observables in the non-resonant domain 1.2 GeV < Mi+i- < Mj/^, taking into 
account experimental bounds on rare B decays such as b —>■ 57. We summarize and conclude 
in Sec. 0. The analytic formulae describing the short-distance effects in the presence of 
SUSY as well as the explicit expressions for the form factors are relegated to the Appendices. 
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II. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL 



In the MSSM, there are new sources of CP violation. In general, a large number of CP- 
violating phases appear in the mass matrices as well as the couplings. After an appropriate 
redefinition of fields one ends up with at least two new CP-violating phases, besides the 
phase of the CKM matrix and the QCD vacuum angle, which cannot be rotated away. For 
instance, in the MSSM with universal boundary conditions at some high scale only two new 
physical phases arise; namely ip^ associated with the Higgsino mass parameter \x in the 
superpotential and ifA connected with the soft SUSY-breaking trilinear mass terms. 

In order to fulfil the severe constraints on the electric dipole moments (EDM's) of electron 
and neutron, one generally assumes that the new phases are less than O(10~ 2 ). Since there 
is no underlying symmetry which would force the phases to be small, this requires fine- 
tuning. Of course, one can relax the tight constraint on these phases by having masses of 
the superpartners in the TeV region; this heavy SUSY spectrum may, however, lead to an 
unacceptably large contribution to the cosmological relic density. 

It has recently been pointed out by several authors that it is possible to evade the EDM 
constraints so that phases of 0(1) still remain consistent with the current experimental upper 
limits. Methods that have been advocated to suppress the EDM's include cancellations 
among different SUSY contributions 031, and nearly degenerate heavy sfermions for the 



first two generations while being consistent with naturalness bounds. The latter can be 



realized within the context of so-called 'effective SUSY' models [25|, thereby solving the 
SUSY FCNC and CP problems. 

To get an idea how supersymmetry affects CP observables in rare B decays, we will 
consider as illustrative examples the following types of SUSY models: 

• MSSM coupled to N — 1 supergravity with a universal SUSY-breaking sector at the 
grand unification scale. 

• Effective SUSY with near degeneracy of the heavy first and second generation sfermi- 
ons. 

In the present analysis, we restrict ourselves to the discussion of flavour-diagonal sfermion 
mass matrices - that is, we assume the CKM matrix to be the only source of flavour mixing.^ 

A. Mixing matrices and new CP-violating phases 

This subsection concerns the mass and mixing matrices relevant to our analysis. In what 



follows, we will adopt the conventions of Ref. 2c 



x One should keep in mind that renormalization group effects induce flavour off-diagonal entries 
in the sfermion mass matrices at the weak scale (see below). 
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1. Charged Higgs-boson mass matrix 



The mass-squared matrix of the charged Higgs bosons reads 

2 _( -B//tan/3 + M^sin 2 /3 + h/vi Bfi + sin (3 cos (3 \ , , 

H± ~ V 5/i + Mjy sin /3 cos /3 By, cot /3 + M ^ cos 2 /3 + t 2 /v 2 ) ' 1 ' 

with 

S/i = ^ sm2P(m 2 Hi + m\ 2 + 2|/i| 2 ), (2.2) 

,o 1 n m?r sin 2 /3 — m?r cos 2 /? 

H 2 = -^| + -^ ^— P -. 2.3 

2 cos 2p 

Here -B and /x refer to the complex soft SUSY-breaking and Higgsino mass parameters 
respectively, m 2 H are the soft SUSY breaking Higgs-boson masses at the electroweak scale, 
and ti j2 stand for the renormalized tadpoles. The mixing angle (3 is defined as usual by 
tan/5 = V2/V1, with t>i j2 denoting the tree level vacuum expectation values (VEV's) of the 
two neutral Higgs fields. In Eq. (|2.2j ) we have adjusted the phase of the \x parameter in such 
a way that Bfi is real at tree level, thereby ensuring that the VEV's of the two Higgs fields 
are real. Consequently, the mass matrix becomes real and can be reduced to a diagonal form 
through a biorthogonal transformation (M^± g ) 2 = OM^ ± T . At the tree level, i.e. £j = 
in Eq. ( |2.1| ), we have 

0=(~ COS f Shl i). (2.4) 
y snip cosp J 

Before proceeding, we should mention that radiative corrections to the Higgs potential induce 
complex VEV's. As a matter of fact, CP violation in the Higgs sector leads to an additional 
phase which, in the presence of chargino and neutralino contributions, cannot be rotated 



away by reparametrization of fields [pTf] . As a result, the radiatively induced phase modifies 
the squark, chargino, and neutralino mass matrices. In the present analysis, we set this 
phase equal to zero. 

2. Squark mass matrices 

We now turn to the 6x6 squark mass-squared matrix which can be written as 

r2 ( Ml r M£„e-^~ 



m !=(m? 9 >* m 2 ' « = U > D > (2 - 5) 

in the (qi, (Jr) basis, and can be diagonalized by a unitary matrix Rg such that 

(M? ag ) 2 = i^M 2 4 ( 2 - 6 ) 
For subsequent discussion it is useful to define the 6x3 matrices 
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(T^) ai = (R q ) ai , (T^ai = (Rg) a , i+3 , q = U,D, (2.7) 

with U and D denoting up- and down-type quarks respectively. Working in the so-called 
'super-CKM' basis || in which the 3x3 quark mass matrices Mjj and Mo are real and 
diagonal, the submatrices in Eq. (|2.5|) take the form 



Mfj = (M?) LL + M 2 + ^M|cos2/3(3-4sin 2 ^)]l, (2.8a) 
lr = Mu\Au - fi* cot (31\, (2.8b) 



Ml = {Ml) RR + M 2 + \m\ cos 2/3 sin 2 9 w l, (2.8c) 



ipfj = &rg(Au - //cot/31), (2.8d) 
M l LL = Vi KM (M 2 ) LL V CK M + M 2 D - cos 2/5(3 - 2 sin 2 B w )% (2.9a) 
M? m = M D |A D -//tan/31|, (2.9b) 

M? m = (M? ) RR + M 2 D - X -M\ cos 2/3 sin 2 d w % (2.9c) 

^ = axg{A D -//*tan/31). (2.9d) 

Here 6*y^ denotes the Weinberg angle, IL represents a 3 x 3 unit matrix, (M~)ll and (M~) RR 
are Hermitian scalar soft mass matrices, and Vckm is the usual CKM matrix. In deriving 
Eq. (|2.9aj ), we have used the relation (M^) LL = V^^M?) ll V C km, which is due to SU(2) 
gauge invariance. Since we ignore flavour-mixing effects among squarks, (M|) n and (M~) RR 
in Eqs. (|2.8|) and ( |2.9|) are diagonal - and hence real - whereas the Aj's are given by 

Ajj = diag(A,, A CJ A t ), A D = diag(A d , A s , A b ), A { = \Ai\e iipA *. (2.10) 



Consequently, the squark mass-squared matrix, Eq. (|2.5|), in the up-squark sector decom- 
poses into a series of 2 x 2 matrices. As far as the scalar top quark is concerned, we have 

2 / m 2 fL + m 2 + \M\ cos 2/3(3 - 4 sin 2 B w ) m t \A t - fx* cot /?|e"^ \ 

* V m t \A t - /I* cot /3\e^i m\ R + m\ + §Mf cos 2/3 sin 2 W J ' 1 j 

where m 2 ^ and m 2 ^ are diagonal elements of {M^) LL and [M^) RR respectively, while (pi can 
be readily inferred from Eq. (|2.8d ). Diagonalization of the stop mass-squared matrix then 
leads to the physical mass eigenstates t\ and t 2 , namely 



tA = f cos e- t sin 9 ie -m \ft L \ = f rg rg« \ / t L 
h J V - sin cos 6 t J V *b / V rgf i«f ; v t R 

where the mixing angle 9$ is given by the expression (—tt/2 ^ di ^ tt/2) 

2m(|A t — /i* cot (3\ 

th tR 



(2.12) 



tan2 ^" W -m? ) + iM 2 cos2/3(3-8sin 2 ^)- (2A3) 



5 



3. Chargino mass matrix 



The chargino mass matrix can be written as 



cos (3 |/i|e* 

where we have adopted a phase convention in which the mass term of the W-ino field, M 2 , 
is real and positive. Note that without loss of generality, we can always perform a global 
rotation to remove one of the three phases from the gaugino mass parameters Mj (i = 1, 2, 3). 

The mass matrix can be cast in diagonal form by means of a biunitary transformation, 
namely 

M dia g = tj* m ^ ±v ^ (2.15) 

where M|± g is diagonal with positive eigenvalues, and U, V are unitary matrices. Solving 
the eigenvalue problem 



( M diag)2 = u*M^Ml ± U T = VM\ ± M^ ± V\ (2.16) 



we find 



/ cos ft/ sin^e-^\ 

1 ~~ y-smOue^ cos 6u )' K ] 

f cosOye-^ 1 sin^e^ (w+< ^ l) \ , * 

\-sm9 v e^v-<t>z) cos9 v e-^ )' [ ' 

with the mixing angles 

tan 2^ = 2 ^ ,mil 2 ^'\; > 1 'r 1 ,: ,; r ; . (2.19) 



2 v / 2M H ,[M|cos 2 /3 


■f \n\ 2 sm 2 (3 + 


|/i|M 2 sin 2/5 cos y^] 1/2 


Mf 


-|/i| 2 -2M 2 / 


cos 2(3 


2 y/2M w [Mf sin 2 (3 - 


f 2 cos 2 (3 + 


\n\M 2 sin 2/3 cos <^] 1/2 


Mf 


-|/x| 2 + 2M 2 / 


cos 2/5 



tan 2ft, = - - - ' 1 — r — , (2.20) 



|/i|sin^sin/5 

tan<^ = - — p , , 2.21 

M 2 cos p + |//| sin p cos 



lul sin cos /3 .„ 

tan^ = - , 7 ^ , (2.22) 

M 2 sin p + cos p cos 



, , sin yyM 2 , sin 2/j 

tan</>1 ~ M 2 (m 2 ± - |/i| 2 ) + |/x|M 2 / sin2/5cos^' 1 } 

Xi 



6 



tan< 



Sm(pJm1 ± — M, 

X.2 



M 2 M^ sin 2f3+\n\ (m 2 ± - Mf ) cos ip^ 

"X-2 



(2.24) 



Here we have chosen — 7r/2 $C 9^ ^ 7r/2, — 7r ^ 0; ^ tt, where i = U,V, and the chargino 
mass eigenvalues read 



m1± = - 

x*2 2 



M| + |/i| 2 + 2M 2 , 



=F {{Ml - |/j| 2 ) 2 + 4M^ cos 2 2(3 + 4M^ [M 2 + | /x | 2 + 2|//|M 2 sin 2/? cos <^]} 



1/2 



^2.25) 



B. SUSY particles and FCNC interactions 

We present here the SUSY Lagrangian relevant to the FCNC processes of interest which 
will also serve as a means of fixing our notation. The interactions of charged Higgs bosons, 
charginos, neutralinos, and gluinos in the presence of new CP phases can be written as [fTTpHfl 



C 



SUSY 
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V2M W 

2 



[cot (3(uMuV CKM P L d) + tanp(uV CKM M D P R d)}H^ 



+ Yl X][u ] {X^P L + Xf R P R )d + u\xf L P L + Xf R P R )l] 

3=1 
4 

+ YA\d\Zk L PL + Z°«P R )d + l\Z^P L + Zi*P R )1\ 
k=i 

8 

- V2g s J2¥ b d\G DL P L - G D *P R )T b d + H.c, 



(2.26) 



where generation indices have been suppressed, and Pl.r = (iTTs)^- The mixing matrices 
in the super-CKM basis are given by 



Xf> = 9 



-VXT Ul + V* 2 T Ur 



Mi 



u 



\/2M w sin/3. 



CKM, 



(2.27a) 



XV» = gU j2 T u ^V CKM 



M, 



D 



V^Mwcos^ 



(2.27b) 



Xf* = -gVZRe, X}* 



gU j2 R u 



V2M w cosp'' 



(2.27c) 



7D L 
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V2 



-N* k2 + \ tan 9 w N* kl ) Y Dl + N* k3 T DR - 
k2 3 kL J ki M w cos(3 



(2.27d) 
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7 D R _ 9 

Zk ""7! 



I tan 6 w N kl T D « + N k3 T D ^ — 

6 Mw cos p 



(2.27e) 



Zk "7! 



(iV^ + tan^iV^r^-iV^r^ 



A2V cos j3 



(2.27f) 





v/2 



2tanfl W 'A r fcl r L « + Y fc3 r Li 



Mr, 



Mw cos /? 



(2-27g) 



(2.27h) 



<^3 being the phase of the gluino mass term M 3 . (For scalar lepton as well as neutralino mass 
and mixing matrices, we refer the reader to Ref. ||26|| .) In the remainder of this section, we 
briefly discuss two SUSY models with quite distinct scenarios for the CP-violating phases. 



C. Different scenarios for the SUSY parameters 



1. Constrained MSSM 



In order to solve the FCNC problem in the MSSM and to further reduce the number of 
unknown parameters, the MSSM is generally embedded in a grand unified theory (GUT). 
This leads to the minimal supergravity (mSUGRA) inspired model, commonly referred to 
as the constrained MSSM 0. In this model one assumes universality of the soft terms at 
some high scale, which we take to be the scale of gauge coupling unification, Mqut, implying 
that (i) all gaugino mass parameters are equal to a common mass M\ji\ (ii) all the scalar 
mass parameters share a common value m ; and (iii) all the soft trilinear couplings are equal 
to Aq. As a result, the mSUGRA model has only two new independent phases which are 
associated with the /i and A parameters. After all, we have at the GUT scale 



tan/5, Mi/2, m , \A \, |/x |, ^a , <Pw> 



(2.28) 



with Mi/2 and mo being real. The parameters at the electroweak scale are then obtained by 
solving the renormalization group equations (RGE's). 

A few remarks are in order here. First, the phases of the gaugino mass terms Mj are 
not affected by the renormalization group evolution, and therefore the low energy gaugino 
mass parameters are real. Second, the phase that appears together with the /x parameter 
does not run at one- loop level so that ip^ = <£> Mo . Moreover, to satisfy the constraints on 
the EDM's of electron and neutron, (p^ has to be of O(10~ 2 ) unless strong cancellations 
between different contributions occur. Third, solving the RGE's for the evolution of the 
CP-violating phase of the A t -term yields a small value for (p^ t , and thus is not constrained 
by the EDM's fl3|,p8]. Lastly, off-diagonal entries occur in the squark mass matrices due to 
renormalization group evolution of the parameters even in the absence of flavour mixing at 
the GUT scale. However, these effects are found to be small and therefore the squark mass 
matrix is essentially flavour diagonal at the electroweak scale (see also Refs. HLT|,p9[|). 
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2. Effective supersymmetry 



As an example of SUSY models with large CP phases, we consider the effective super- 
symmetry picture |ES| without assuming universality of sfermion masses at a high scale. 



Within such a framework, the first and second generation sfermions are almost degenerate 
and have masses above the TeV scale, while third generation sfermions can be light enough 
to be accessible at future hadron colliders. Consequently, FCNC reactions as well as one-loop 
contributions to the EDM's of electron and neutron are well below the current experimental 
bounds. 

However, it should be noted that the EDM's also receive contributions at two-loop level 
involving scalar bottom and top quarks that may become important for phases of order 
unity in the large tan/3 regime ||30|| . In our numerical work, tan/3 is assumed to be in the 
interval 2 ^ tan (3 ^ 30. 

III. RARE B DECAYS AND NEW PHYSICS 

A. Short-distance matrix element 

Let us start with the QCD-corrected matrix element describing the short-distance inter- 
actions in b — > s(d)l + l~ within the SM. It is given by 



Gfoi. 



2cf 



cf-c 10 )(H(k)\fj,P L b\B(p)) 



q 2 



(H(k)\fia^q» (m b P R + m f P L ) b\B(p)) 



I^PlI + (cio - -c 10 )T-fP R l }, (3.1) 



where q is the four-momentum of the lepton pair, and H = K, K* (it, p) in the case of / = s 
(f = d). In the SM, the Wilson coefficients cf and c±o are real with values of —0.314 and 
—4.582 respectively, and the leading term in Cg ff has the form [ |3"T| , |3"2] 



cf = c 9 + (3ci + c 2 ) jg(m c ,g 2 ) + X u g(m c ,q 2 ) - g(m u ,q 2 ) |h , (3.2) 

where eg = 4.216. The Wilson coefficients will be discussed in detail in Appendix [B]. In the 
above expression 

v v* { - x2 (p- ir i) for / = s > 

K = « Pi} ~P)-V 2 _ iV f f _ . (3.3) 

^ [ (l_ p )2 + r? 2 (1 _ p )2 +r? 2 IOT / 

with p and r\ being the Wolfenstein parameters |3"3]|, where the latter reflects the presence 



of CP violation in the SM. For definiteness, we will assume p = 0.19 and rj = 0.35 [34| 
Finally, the one- loop function g{m^ q 2 ) at the scale p R = m b is given by0 



2 In order to avoid confusion with the fi parameter of the superpotential, we use the notation /ir 
for the renormalization scale. 
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g 8 4 2 

g(m h q 2 ) = -- hx{mi/m h ) + — + -y t - ~(2 + ^)a/Ti - & 



x 6(1-1/,) 



1-VT^ 



0(1/, - 1)2 arctan ^ , (3.4) 



where = 4mf/q 2 . Observe that cc and mm loops provide absorptive parts that are manda- 
tory, as we show below, for a non-zero partial width asymmetry besides the presence of a 
CP-violating phase. 

B. Wilson coefficients and new physics 

Throughout this paper, we will assume that in the presence of new physics there are 
no new operators beyond those that correspond to the Wilson coefficients appearing in 
Eq. ( |3.1|) . (For a discussion of the implications of new operators for rare B decays, see, e.g., 
Ref. ||35|| .) Thus, the effect of new physics is simply to modify the matching conditions of 
the Wilson coefficients, i.e. their absolute values and phases at the electroweak scale. 

As a result, we are left with additional SUSY contributions at one- loop level to the 
Wilson coefficients cf , Cg ff , and cio in Eq. Q3.1J) . In fact, they arise from penguin and box 
diagrams with (i) charged Higgs boson up-type quark loops; (ii) chargino up-type squark 
loops; (iii) neutralino down-type squark loops; and (iv) gluino down-type squark loops. 
Thus, the short-distance coefficients can be conveniently written as 

Cl (M w ) = cf M (M w ) + c?\M w ) + cf{M w ) + cf{M w ) + 4(M W ) (i = 7, . . . , 10). (3.5) 

The explicit expressions for the various Wilson coefficients are given in Appendix |B[ Since 
we limit our attention to flavour-conserving effects in the squark sector, the neutralino and 
gluino exchange contributions in Eq. ( |3.5| ) will be omitted in our numerical calculations. 
For future reference, we parametrize the new physics contributions as follows: 

Ci (M w ) = , R ,* i ^8-1 , } 

where x is rea l to a good approximation within the models under study. 



C. Resonant intermediate states 



We have considered so far only the short-distance interactions. A more complete ana- 
lysis, however, has also to take into account resonance contributions due to mm, dd, and 
cc intermediate states, i.e. p, u, J/ip,ip', and so forth. A detailed discussion of the various 
theoretical suggestions of how to describe these effects is given in Ref. p5 |. 



We employ here the approach proposed in Ref. |57| which makes use of the renormalized 



photon vacuum polarization, Il^ ad (s), related to cross-section data[] 



This method assumes quark-hadron duality and rests on the factorization assumption. 
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„ ( * _ ^tot(e + e -> hadrons) 

^had(s) = r— T _> , 3.7 

cr(e + e — > ) 

with s = q 2 . In fact, the absorptive part of the vacuum polarization is given by 

lmlP bad (s) = ^R bBd (s), (3.8) 



whereas the dispersive part may be obtained via a once-subtracted dispersion relation [38 



RelTUs) = / (3.9) 



4A/2 



with P denoting the principal value. For example, in the case of the J/ip family 
(i.e. J/ip, ip', . . .) the imaginary part of the one-loop function g(m c ,s), Eq. ( |3.4| ), can be 
expressed as 



lmg(m c , s) = ee it£ nt (s) + R^(s), (3.10) 

where the subscripts 'cont' and 'res' stand for continuum and resonance contributions re- 
spectively, while the real part is given by 

8 4 s °T R^^fs') 

Reg(m c ,s) = --Hm c /m b )-- + -P / h f 1 U s'. (3.11) 

9 9 6 J s'{s' — s) 

The contributions from the continuum can be determined by means of experimental data 



given in Ref. |39|| , whereas the narrow resonances are well described by a relativistic Breit- 
Wigner distribution. 

However, in order to reproduce correctly the branching ratio for direct J/ip production 
via the relation (V c5 = J/ip,ip', . . .) 

B(B -> HVcb -> Hl + r) = B(B -> HV c5 )B(V c5 -> l+l~), (3.12) 

where H stands for pseudoscalar and scalar mesons, one has to multiply R/ls in Eqs. (|3.10 ) 
and ( p.ll| ) by a phenomenological factor k, regardless of which method one uses for the 
description of the resonances [40, 41] .0 Using the form factors of Ref. [43] (see next section) 
together with experimental data on B(B — > K^J/ip), B(B —> K^ip'), and B(B~ — > 
7r~J/ip) |44]], we find a magnitude for k of 1.7 to 3.3. At this point two remarks are in order. 
First, the branching ratio for direct J/ip and ip ; production is enhanced by a factor k 2 , while 
it is essentially independent of k outside the resonance region. Second, the numerical results 



4 Strictly speaking, it is the term (3ci + C2)RrL - in the approximation of Eq. 
to be corrected to k(3ci + C2)Rrls , taking into account non-factorizable contributions in two-body 
B decays (see, e.g., Ref. J|2j). 



3.2) - which has 
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for average CP asymmetries in the non-resonant continuum 1.2 GeV < ^/s < 2.9 GeV are 
not affected by the uncertainty in k. 

Similar considerations also hold for uu and dd systems except that the p resonance is 
described through 

1 / 4M^ 3/2 



RPJ S ) = -[1-—IL) \F n (s)\ 2 , (3.13) 



where the pion form factor is given by a modified Gounaris-Sakurai formula E5 



IV. THE DECAYS B -> K"W l+l~ AND B -> jr(p) l+l~ 

The hadronic matrix elements in exclusive B decays can be written in terms of q 2 - 
dependent form factors, where q 2 is the invariant mass of the lepton pair. In the work 
described here, we employ heavy-to-light B — ► and B — ► n(p) form factors determined 
by Melikhov and Nikitin within a relativistic quark model. To get an estimate of the 
theoretical uncertainty that is inherent to any model for the form factors, we also utilize the 



parametrization of Colangelo et al. [IB], which makes use of QCD sum rule predictions. 

For simplicity of presentation, we do not display corrections due to a non-zero lepton 
mass, which can be found in Refs. [p2], fl7]l . (The same remark applies to the light quark 
masses m s ,d-) Henceforth we shall denote pseudoscalar and vector mesons by P = K, tt and 
V = K*, p respectively. 



A. B — > P transitions 

1. Form factors 

The hadronic matrix elements for the decays B — > P can be parametrized in terms of 
three Lorentz-invariant form factors (see Appendix [D]for details), namely 

(P{k)\j lfl P L b\B(p)) = h [ {2 V -q)^ + {q 2 )+q,f_{q 2 )}, (4.1a) 

(P(k)\fia^P L!R b\B(p)) = ~[{2p - q),q 2 - (M 2 B - M 2 )q,}s(q 2 ), (4.1b) 

with q = p — k. Here we assume that the form factors are real, in the absence of final-state 
interactions. Note that the terms proportional to q^ may be dropped in the case of massless 
leptons. 



2. Differential decay spectrum and forward-backward asymmetry 
Introducing the shorthand notation 

X(a, b, c) =a 2 + b 2 + c 2 - 2(ab + bc + ac), (4.2) 
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(4.3) 



and recalling s = q 2 , the differential decay rate can be written as (neglecting mi and rrif 
dT(B -> PZ+Z-) G^a 2 



ds d cos Q] 



4 ti f + (s) + 2c?m b s( S )\ 2 + \c 10 f + (s)\< 



sin 2 0,. 
(4.4) 



Here 8i is the angle between l~ and the outgoing hadron in the dilepton centre-of-mass 
system, and the Wilson coefficients are collected in Appendix |B[ Defining the forward- 
backward (FB) asymmetry as 



A FB {s) 



\ Q dT 

dcos9l l~l 7T 

o as a cos Ui 



a cos vi - 

•l 



' lis d cos 9i 



d cos 9i— — 

as a cos 



dr 

dcos6l l~l 7T 

4 as a cos 



(4.5) 



which is equivalent to the energy asymmetry discussed in Ref. Jll[, it follows directly from 
the distribution in Eq. ( [4.4[ ) that Afb vanishes in the case of B — > Pl + l~ transitions. We 
note in passing that, given an extended operator basis (e.g. in models with neutral Higgs- 
boson exchange), new Dirac structures // and Z75Z may occur in Eq. ( |3.1| ), giving rise to a 
non-zero FB asymmetry in B — > Pl + l~ . 



B. B 



V transitions 



1. Form factors 

The hadronic matrix elements describing the decays B — ■> V are characterized by seven 
independent form factors, which we present in Appendix [D], defined through (eoi23 = + 1) 

(V(k)\h,P L b\B(p)} 

= ie^e v *p a qPg(q 2 ) - \^J{q 2 ) + (e* • q)[(2p - q),a + (q 2 ) + g M a_(g 2 )]}, (4.6a) 



(V(k)\fia, v q»P R , L b\B(p)) = ie, ua(3 e»*p a qPg + (q 2 ) T ^*[g + (q 2 )(M 2 B - My) + q 2 g_{q 



1 



±-{e*-q)\{2p-q) li 



g + (q 2 ) + -q 2 h(q 2 ) 



9-{q 2 )-^M 2 B -M 2 )h{q 2 ) 



(4.6b) 



being the polarization vector of the final-state meson, and q = p — k. 
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2. Differential decay spectrum and forward-backward asymmetry 



The differential decay rate for B — > Vl + l in the case of massless leptons and light 
quarks takes the form (/ = s or d) 



dT{B -> W+Z- 



GW 



2 9 tt 5 M| 1 46 4/1 



|F t6 F t * f A(s) + B(s) cosOi + C(s) cos' 6i 



ds d cos 6i 

with Xy as in Eq. (13). The quantities A, B, and C are 



A(s) 



2X1 

Ml 



1 / 2sM 2 

sM v a x {s) + - (^1 + )«2(s) +4a 3 (s) ■: (/> • ( /)n ](.s) 



(4.7) 



(4.8) 



S(s) = 8X v Re{ct [cf sA^ - c?m b (A x B y + A y B x )}} 



C(s) 



2X1 

M v 



sMy«i(s) - ^a 2 (s) - X v a 3 (s) - (fc ■ q)a A (s) 



where k ■ q = (M| - M$ - s)/2, and 
ai(s) = (|cf| 2 + |c 10 | 2 )^ + 



2 , , 2 , 4j cf | 2 m 2 p2 4Re (cfcf*)m fc 

Z -D,, 



^4 a- B x , 



(4.9) 
(4.10) 



(4.11a) 



a 2 {s) = ai(s) s _^, a 3 (s) = <*i(s) x - 



(4.11b) 



a 4 (s) = (|cf | 2 + |c 10 | 2 )AA + Uli^S^, - 2Reic \ cV)m \ A y B z + A z B y ), (4.11c) 



in which the A^s and B^s are defined as 

A x = g(s), A y = f(s), A z = a+(s) 



(4.12a) 



B x = g + (s), B y = g + {s){M 2 B - My) + sg^s), B z 



9+( s ) + 2 s/l ( s ) 



(4.12b) 



The complex Wilson coefficients cf , Cg ff , and c\q are given in Appendix |B[ 

Finally, using Eqs. ( f4.5|) and (|4.7| ), we derive the forward-backward asymmetry 



ReKJcgsA^ - cfm b (A x B y + A y B x )]} 
[3A(s) + C(s)} 



(4.13) 
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C. CP- violating observables 



To discuss CP- violating asymmetries, let us first recall the necessary ingredients. Suppose 
the decay amplitude for B — > F has the general form 

A(B -»• F) = e^Axe^ 1 + e^ 2 A 2 e iS \ (4.14) 

where 5i and 4>% denote strong phases (CP-conserving) and weak phases (CP- violating) re- 
spectively (Ai and A 2 being real). Together with the decay amplitude for the conjugate 
process 

A(B -»• F) = e-^Axe** 1 + e"^ 2 A 2 e iS \ (4.15) 

which can be obtained from Eq. (|4.14|) by means of CPT invariance, we may define the CP 
asymmetry as 

|^l| 2 -|g| 2 = -2rsin0sin£ 
CP |.A| 2 +|.A| 2 1 + 2rcos0cos<5 + r 2 ' l ' J 

with r = A2/A1, = 0i — 02, and 5 = 5i — 5 2 . As can be easily seen from Eq. (|4.16|) , a 
non-zero partial rate asymmetry requires the simultaneous presence of a CP- violating phase 
as well as a CP-conserving dynamical phase 5, the latter being provided by the one-loop 
functions g(m c ,s) and g(m u ,s) present in the Wilson coefficient cf 1 [Eq. (|3.2|) 1. Notice that 
in the limit in which the charm quark mass equals the up quark mass there is no CP violation 
in the SM. 

Given the differential decay distribution in the variables s and cos 9i, we can construct 
the following CP-violating observables: 



d cos 9, 



dT d iff 

1 Ac. A ~ZZ a. r rl rO 



A& { .) = # "^L, I = / =f / . (4.17a) 



where we have introduced 



diflf 



dcose l ry m \ JD ' S J0 J ~ l 

s dsd cos 9i 



T(B -> Hl + r) - T(B -> Hl + r), (4.17b) 



r sum = r(B -> m + r) + r(s -> Fz+r), (4.i7c) 

with H = K(*\ it, p. It should be noted that the asymmetry A^ p represents a CP- violating 
effect in the angular distribution of l~ in B and B decays while Aq P is the asymmetry in 
the partial widths of these decays. As can be seen from Eqs. (|4.4|) , (|4.7|) , and (|4.13|) , the 



latter involves the phases of and Cg ff while the former is also sensitive to the phase of the 
Wilson coefficient cio- 

V. NUMERICAL ANALYSIS 

Given the SUSY contributions presented in the preceding sections, we now proceed to 
study the implications of supersymmetry for exclusive B decays. 
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A. Experimental constraints 



In our numerical analysis, we scan the SUSY parameter space as given in Ref. |D| and 
take as input the parameters displayed in Appendix In addition, we take into account 
the following experimental constraints: 

• From the measurement of the inclusive branching ratio B(B — ► Xgj), which probes 



c% \, one can derive upper and lower limits 48 



2.0 x 10" 4 < B{B -> X sl ) < 4.5 x 1(T 4 (95% C.L.; 



(5.1) 



This is specially useful to constrain extensions of the SM. Indeed, following the model- 
independent analysis performed in Ref. |[49|| , and taking the Wilson coefficient c| in 
leading- log approximation [see Eq. flB2"| ) of the Appendix], we obtain 



B{B -> X sl ) « [0.801 + 0.444|i2 7 | 2 + 0.002|i? 8 | 2 + 1.192R ej R 7 
+ 0.083Re J R 8 + 0.061Re( J R 7 i^)] x 10" 4 , 



(5.2) 



where _R 7 and Rs are as in Eq. (|3.6|) . From Fig. |l] we infer that the present b — > 57 
measurement already excludes many solutions for R7. 

A CDF search for the exclusive decays of interest yields the upper limits B(B° — > 
K*°n + fi-) < 4.0 x 10~ 6 and 23(5+ -> K+fi+fi') < 5.2 x 10" 6 at the 90% C.L. |U[. 
Note that the K*°p: + pi~ upper limit is close to the SM prediction ]19, 51]. As for the 
modes B — > 7r/ + Z~ and 5 — > pl + l~, we are not aware of any such limits. 

Non-observation of any SUSY signals at LEP 2 and the Tevatron imposes the following 



lower bounds [44] : 



m^± > 86 GeV, m s > 43 GeV, m h > 86 GeV, m g - > 260 GeV, m H ± > 90 GeV. (5.3) 



B. CP asymmetries 

As mentioned earlier, we investigate CP asymmetries in the low dilepton invariant mass 
region, i.e. 1.2 GeV < yjs < (Mj/^p — 200 MeV), which is of particular interest because 
the low-s region is sensitive to the Wilson coefficient Cj (in the case of B — > Vl + l~). In 
fact, it can receive large SUSY contributions and be complex, as well as change sign, while 
being consistent with the experimental measurement of b — > 57. On the other hand, the 
new-physics effects are known to alter the remaining Wilson coefficients Cg ff and c±q only 
slightly within mSUGRA and effective SUSY with no additional flavour structure beyond 
the usual CKM mechanism |5,12,19]. Moreover, CP asymmetries above the J /if) resonance 



are dominated by cc resonant intermediate states, whereas below 1 GeV the p resonance has 
a strong influence on the asymmetry. This can be seen from Fig. |^, where we show the CP 
asymmetry [Eq. (|4.17a )] between B~ — > tt _ / + /~ and B + — > tt + Z + /~ as a function of 



the dilepton invariant mass within the SM and in the presence of SUSY contributions with 
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123456 123456 

07 07 

FIG. 1. Allowed region for |i?7| and the corresponding CP-violating phase <j)j as determined 
from the inclusive measurement of b — ► S7 rate, using the leading-log expression for Cy ff . Diagrams 
(a) and (b) correspond to different values of x [Eq. (|3.6|)]. 



new CP-violating phases. It is evident that the predictions for CP asymmetries suffer from 
large theoretical uncertainties in the neighbourhood of the p resonance and above the J/i/j, 
as discussed in Sec. |T|. 

Using Eqs. ( |4.4| ) and ( |4.7| ) together with the definition for CP-violating asymmetries, 
Eqs. ( |4.17| ), we can summarize our main findings as follows: 



1. CP violation in B — > Pl + l 

• The CP-violating asymmetry in the l~ spectra of B and B decays, A^-p, vanishes in 
the case of B — > P transitions. Within the framework of the constrained MSSM with 
phases of O(10~ 2 ) numerical values for the average asymmetry (^4cp) m the low-s 
region are comparable to the SM predictions with asymmetries of 0.1% and —5% for 
b — > s and b — > d transitions respectively. 

• Our results for A^ p between the decays B — > Pl + l~ and B —>■ Pl + l~ in the context of 
effective SUSY with a light stop ti and phases of Oil) are shown in Figs. || and | for low 
and large tan/? solutions which correspond to Rei?7 > and KeR 7 < respectively. 
Observe that the CP-violating asymmetry Aq P in B — >• P depends only weakly on the 
sign and phase of Cy ff . This is due to the fact that Cy ff , which is constrained by the 
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FIG. 2. CP-violating partial width asymmetry ^4^p m the decays B~ — » ir~l + l~ and 
B + — > vs vi, s = s/Mg, including p, u, and J /if), if/, etc. resonances, and employing 

the form factors of Ref. E3|. (a) Within the SM and (b) in the presence of new CP phases and 
a real CKM matrix. For the sake of illustration, we have chosen \Rj\ = 1.6, |i?gio| = 1, X = 1-5, 
07 = 7r/2, and 4>9,io = 0.01. 



b — > 57 measurement and not enhanced by a factor of 1/s in the low-s region, is nearly 
one order of magnitude smaller than the leading term in Cg ff [cf. Eq. (3.2)]. 



The CP asymmetry in the partial widths of B — > Kl + l~ and B — > Kl + l~ changes 
sign for large values of 0g, while \A<L V \ ^ 1% (see Fig. |3|). However, non-standard 
contributions to <pg are found to be small and hence Aq F ~ O(10~ 3 ). On the other 
hand, average asymmetries of —(5—6)% are predicted for (A^p) ^ n ^ ne case of B ^ n, 
even for values of <pg as small as 10~ 2 (see Fig. |]). Given a typical branching ratio of 
10 -8 and a nominal asymmetry of 6%, a measurement at 3a level requires 3 x 10 11 bb 
pairs. (This rather challenging task might be feasible at LHC and the Tevatron.) 

The small magnitude of the CP asymmetry is also due to a suppression factor mul- 
tiplying the indispensable absorptive part in Cg ff , which is only slightly affected by 
new-physics contributions. Indeed, it follows from Eqs. ( |3.2j) and (|4.16|) that 

(3d + C2) 

Aqp oc r sin (ft sin 5 ~ : — : sin <p sin 8 ~ 1CT 2 sin <\> sin 5, (5.4) 

M 

where the weak and strong phases <ft and S can be of order unity. 
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3.6 



FIG. 3. CP-violating partial width asymmetry A^ p in the decays B — ► Kl + l and B — > 
as a function of 0g and |i?7| for a dilepton invariant mass of s = 4GeV 2 within effective SUSY, 
(a) tan/3 = 2 with 4> 7 = 0.4, |i? 9 | = 0.96, \R 10 \ = 0.8. (b) tan/3 = 30 with 7 = 2.5, |i? 9 | = 0.99, 
\Rio\ =0.9. 
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3.6 

FIG. 4. CP-violating partial width asymmetry ^4^p between B~ — > ir~l + l~ and i? + — > 7r + Z + Z — 
as a function of <pg and |i?7| for s = 4GeV 2 within effective SUSY, (a) tan/3 = 2 with <f>j = 0.4, 
\R 9 \ = 0.96, |i?i | = 0.8. (b) tan/3 = 30 with (f> 7 = 2.5, |i? 9 | = 0.99, |i?i | = 0.9. 
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Numerical estimates for average CP asymmetries are mildly affected by the parame- 
trization of form factors (see also Refs. [22, 52||). 



2. CP violation in B — > Vl + l~ 

• The contribution of the Wilson coefficient Cy ff (or equivalently -R7) to the decay rate 
in B — > V modes is enhanced by a factor of 1/s in the low-s region. As seen from 
Fig. |5], in the case of B — > p, the CP-violating asymmetry Aq F can change sign for 
tan/3 = 2 (i.e. Rei?7 > 0), while for large tan/5 it is always negative. For small values 
of 4>g an average CP asymmetry of about —5% is predicted for both tan (3 = 2 and 
tan (3 = 30 solutions. Since the distributions of Aq P for B — > K* are very similar to 
the ones obtained for B — ► K, we refrain from showing the corresponding plots. 

• As we have already mentioned, the CP-violating asymmetry in the angular distribution 
of l~ in B and B decays can, in principle, probe the phase of the Wilson coefficient 
cio- This is shown in Fig. |6], where we have plotted the CP asymmetry as a function of 
09 and 0io for large tan/3 (i.e. Rei? 7 < 0). Unfortunately, the mSUGRA and effective 
SUSY predictions for the average asymmetry (Aqp) t urn ou ^ to be unobservably small. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the consequences of new CP-violating phases for ex- 
clusive B decays within the framework of supersymmetric extensions of the SM, ignoring 
inter generational mixing in the squark sector. We have examined CP- violating asymmetries 
in the partial widths as well as angular distribution of l~ between the exclusive channels 
b — > s(d)l + l~ and b — > s{d)l + l~ in the invariant mass region 1.2 GeV < M1+1- < 2.9 GeV. 
The essential conclusion of our analysis is that it is not sufficient to have additional CP 
phases of 0(1) in order to obtain large CP- violating effects. 

Within the constrained MSSM and effective SUSY with a complex CKM matrix and 
additional CP phases, we obtain values for the average asymmetry (Aq P ^ of about —6% 
(—5%) in the decays B — > n(p)l + l~ and B — > 7f(p)/ + /~, taking into account experimental 
constraints on EDM's of electron and neutron, as well as rare B decays such as b — > 37. 
As for the asymmetry in the angular distribution, (A^p^, it probes the phase of the Wilson 
coefficient Cio, but will be unobservable at future colliders. Numerical estimates of the CP 
asymmetries in the decays B — > K^*H + l~ and B — > R(*H + l~ turn out to be small (less than 
1%) and are comparable to the SM result. 

Our analysis shows that the smallness in CP asymmetries is mainly due to the coefficient 
(3ci + c 2 )/|c 9 | which multiplies the requisite absorptive part in Cg ff [Eq. (|3.2j) 1, and which is 
only slightly affected by the new-physics contributions discussed in Sec. [TIT]. Therefore, any 
sizable CP-violating effect in the low-s region requires large non-standard contributions to 
the short- distance coefficient Cg ff and/or cio, as well as additional CP phases of 0(1). By the 
same token, any large CP-violating effect would provide a clue to physics beyond the SM. 
A detailed discussion of this point will be given elsewhere. 

One could argue, however, that the inclusion of flavour off-diagonal contributions 
(i.e. gluino and neutralino diagrams) to the Wilson coefficients might lead to higher CP 
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3.6 



FIG. 5. CP asymmetry in the decays B — > p l + l and B + — > as a function of 

(j)g and |i?7| for a dilepton invariant mass of s = 4GeV 2 within effective SUSY, (a) tan/3 = 2 with 
07 = 0.4, |i? 9 | = 0.96, \R W \ = 0.8. (b) tan/3 = 30 with (j> 7 = 2.5, |i2g| = 0.99, |i?i | = 0.9. Note 
that (a) and (b) correspond to Re Rj > and Re R-j < respectively. 
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FIG. 6. CP-violating asymmetry A^ p between (a) B —> K*l + l~ and B — * K*l + l~, and (b) 
B~ — > p~l + l~ and B + — * p + l + l~ for large tan/3 as in Fig. ||. 
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asymmetries. In fact, it has been pointed out in Refs. |29|,|53| that even in the presence of 
large supersymmetric CP phases, a non-trivial flavour structure in the soft-breaking terms 
is necessary in order to obtain sizable contributions to CP violation in the K system and to 
CP asymmetries in two-body neutral B decays (see also Ref. p4j). Using the mass insertion 
approximation, such effects have recently been studied in Ref. which predicts a partial 
width asymmetry for b — > sl + l~ of a few per cent in the low-s domain. 

Finally, we would like to recall that, nevertheless, large CP asymmetries may occur in 
rare B decays like the observed b — ► 57 modes, where A C p can be substantial (up to ±45%) 
in some part of the parameter space [T4j, 
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APPENDIX A: NUMERICAL INPUTS 

Unless otherwise specified, we use the experimental values as compiled by the Particle 
Data Group |44j and the parameters displayed in Eq. (|A~1|). 



m t = 175 GeV, m b = 4.8 GeV, m c = 1.4 GeV, m s = 170 MeV, m d = 10 MeV, ( 
m u = 5 MeV, a = 1/129, A QCD = 225 MeV. ^ ' 



APPENDIX B: WILSON COEFFICIENTS AND SUSY 

For the sake of convenience, we provide in this appendix formulae for the Wilson coeffi- 
cients Cy ff , Cg ff , and c w in the presence of SUSY, using the results derived in Refs. |p|, pTT| , |3T] , |55 



Since we study the case of massless leptons, we retain only those contributions that do not 
vanish in the limit mi — > 0. As for r leptons in the final state, there are further charged and 
neutral Higgs-boson contributions [see also Eqs. ( p. 271 )]. 
Introducing the shorthand notation 



_ a s {M w ) ^ rw= "±^ r ^ ± = ^L ; t a = ^a^ (B|) 
a s (m b ) ' My/ m //±' B m V 

the Wilson coefficient evaluated at = rrib has the form (in leading-log approximation) 

cf = vl 6/23 c 7 (M w ) + |(^ 23 - V f)c 8 (M w ) + hit/, (B2) 



i=i 



with the coefficients a i; hi tabulated in Ref. JHj]. Recalling Eqs. ( J2.27[ ) and (|3.5|), and using 
the one-loop functions /j listed in Appendix |C|, the various contributions to cj : % (Mw) can 
be written as follows: 
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Standard model: 



c* M (M w ) = - A r w h{r w ). 



(B3) 



Charged Higgs boson: 



c^(M w ) = —[r H ±h{r H ±) cot 2 (3 + 2f 2 {r H± )}. 



(B4) 



Chargino:[] 



4 (M, 



j 6 2 

EE 



y *6 */ a=l j=l X f 



X 



(^WU(^) -2(x^ t ) i;(; cv; "•) 

where we have defined 



m-± 



n 



1 for / = d, 

2 for / = s. 



(B5) 



(B6) 



For completeness, we also give the expressions for the neutralino and gluino contributions 
which vanish in the limit of flavour- diagonal squark mass matrices. 



Neutralino: 



cf{M w ) 



1 b 4 M 2 

1 \^ \^ M W 

Qg 2 V fb V; f ^ ^ m 2 



• (B7) 



• Gluino: 

Jim \- 4 ft 2 y M w 



y rrib y 



(B8) 



The corresponding expressions cf M (M w ), . . . , c| (M w ) are obtained changing f\ — > gi in 
Eqs. ( |B3| )-(P7|), with gi collected in Appendix [CJ while the gluino contribution reads 



5 Notice that the one- loop function appearing in the last term of Eq. ( P5| ) is actually f 2 + 5/2. 
However, using the explicit form for the squark mixing matrices [Eqs. ( |2.27| )1, the constant term 
vanishes - reflecting the unitarity of the mixing matrices. 
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Till, 



d a ^ 



(B9) 



The Wilson coefficient Cg at = m b in next-to-leading approximation is given by 

a s (m b ) 



Cg — Cg 



1 + 



7T 



-u{s/m b ) 



+ #(m c , s) (3ci + c 2 + 3c 3 + c 4 + 3c 5 + c 6 ) 



+ A u [ff(m C) s) - g(m u , s)} (3ci + c 2 ) - ^g{m f , s) (c 3 + 3c 4 ) 

1 2 
- -</(m 6 , s) (4c 3 + 4c 4 + 3c 5 + c 6 ) + - (3c 3 + c 4 + 3c 5 + c 6 ) 



(BIO) 



where A u and g(rrii,s) are defined in Eqs. ( |3.3| ) and ( ^.4[ ) respectively, with s = q 2 . As far 
as the Wilson coefficients Ci-c 6 are concerned, we have numerically 

ci = -0.249, c 2 = 1.108, c 3 = 0.011, c 4 = -0.026, c 5 = 0.007, c 6 = -0.031, (Bll) 

using the values given in Appendix |A]. Further, 



c 9 = c 9 (M w ) - - + P + P E J2 Ei 



(B12) 



with i = SM, H ± , x ± , X°i 9i an d 



c 9 {M w ) = 



sin 2 6 W 



AZ l + 



(B13) 



where the analytic expressions for P , Pe, and E b are given in Ref. [31]. Since Pg -C P c 



o- 



we shall keep only the SM contribution in the last term of Eq. ( [B12[ ). Moreover, as discussed 
in Ref. ||T| , the order a s correction in Eq. ( piO|) due to one-gluon exchange may be regarded 
as a contribution to the form factors, and hence we set to = in Eq. ( B10|) . 
Turning to the Wilson coefficient Cio, it has the simple form 



cio(M, 



w 



^ sin 2 e 



ir 



(B14) 



Note that the corresponding operator does not renormalize and thus c w (M w ) = Ci (/ii?) 
The expressions for the various contributions to Y and Z read as follows:[] 

• Standard model: 



6 Regarding the expressions for the chargino and neutralino box-diagram contributions, and the 
sign discrepancy between Ref. || and Ref. [11|, we confirm the results of the latter. 
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• Charged Higgs: 



Y SM = \Mr w ), Z SM = ±f 10 (r w ). 



Y H± =Y Z H \ Z H± = zf + zf, 



7H± 



~U{r H± ) cot 2 /?, Y Z H± = Z$ ± = -\r w h{r H± ) cot 2 (3. 



• Chargino: 



(B15) 

(B16a) 
(B16b) 

(B17a) 



1 M 2 X s " 

^t/ a=1 , =1 m u a 



(B17b) 



6 2 



F # - Z Z ~ 9 2 V V * H\ ( X f if W^7 L )b3 
Z9 V tb V tf a ,b=li,j=l < 



1 



- C2(ml a ,m 2 -±,m 2 ~±)5 ab V* 1 V jl + -m^±m^±c (ml a , m^,m\±)b ah U iX U] x 



c 2 (m 2 ±! m 2 a ,m 2 J(r^r^t )afe% 

(B17c) 



b ° X <? 2 ^ 



6 2 



rEE^^-^f^as^K^m^^l^^J^^!], (B17d) 



t& t/ a=l i,j=l 

with mjj^m^) in the case of e + e~(/i + /i~) in the final state. 
• Neutralino: 

Yf^Yf + Y*, Z*° = Z? + zf + Z } 



x 

box' 



(B18a) 



< = -^m, 1 1 ^y^u^uh^i (bib,) 

y *& t/ a =l fc=l d a 



^ = < 



1 6 4 { Dif 

o„2T/ T/* S 1 ( Z k L )na(Z° L ) b 3 

19 v tb v t f „ b=lk , =1 L 



x 



- im^.m^.c (m| a ,m 2 o,m 2 o)(5 afe (Ar fc 3Ar* - iV fe4 ^* 4 ) 



(B18c) 
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box - * bm ^W^box -I- o„2T/ T/* / / ^ U A * MoV^/ 



/ a=l fc,/=l ^ 



x d 2 (m|o, m| ? , m| a , m? 2 )(A^* 2 + tan^^i)(^2 + tan0 w iVu) 

+ ^m^o rriyo d (m|o , m|p , to| q , to? a ) (N k2 + tan W iV** ) (N* 2 + tan W N* x ) 



z & = otE E Uz^U^U^e 

9 v tb v tf tt=1 w=1 L 



d2 ( m l°' m l?' TO !> m i, 5 ) N ki N u + -TO^TO f odo(m|o, to|o, to| q) TO? jAAfci^ 
with m|- i4 (m ; - 25 ) for e + e~(yU + /i~) in the final state. 



Gluino: 



,2 6 11^2 



n 2 M l 
a to tj a= i d a 



Yl = 4 = E [(^% a (^) &3 C2(TOi ) TOl a ,TO|)(^G^t ) 



,,6=1 

The functions /j, q, and are given in Eqs. (|CT|)-(|C4|) below. 

APPENDIX C: AUXILIARY FUNCTIONS 

Here we list the functions /j, g i: c iy and di introduced in the previous section. 

„ . . -7 + 5x + 8x 2 x(2 - Sx) , 
A(X)= 6(1 -x)3 TT^j 4 " ' 

„ . . a;(3 — 5a;) x(2 — 3a;) 



, , . 2 + 5x - x 2 x . 

/3( " )= 6(1 -x)3 + (T^F lnX ' 



/4(x) = 27^ + (T^F lna; ' 
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fb(x) = + . X , 2 lux, 

1 — x (1 — xy 



, , , x(38 - 79a; + 47a; 2 ) x(4 - 6a; + 3a; 3 ) , 
h\ x ) = — ^ + t, —a lnx , 



6(1 -a;) 3 



(1-xY 



r / \ 52- 10 la; + 43a; 2 6 - 9x + 2a; 3 , 
/tW = — — + —T, ^2— lnx ' 



6(1 -a;) 3 



(1-xY 



r . . 2 - 7x+ llx 2 6a; 3 . 

f8iX)= (1-X)3 + (T^)I lna; ' 



g(4-a) 3a: 2 
1 — x (1 — x)^ 



p , , x(108 - 259a; + 163a; 2 -18a; 3 ) 8 - 50a; + 63a; 2 + 6a; 3 - 24a; 4 , 
fmx) = — T , — lux, 



2(1 - a;) 3 



(1 - a;) 4 



gi{x) = -3/ 3 (x), g 2 (x) = 3[/ 4 (x) - 1/2], # 3 (x) = -3/ 3 (x), # 4 (x) = -3/ 4 (x), 



#>(x) = g 



11 -40a;- 19a; 2 3x(l - 9x) n 

- + — 7t rr^lna; 



2(1 -xf 



(1-xY 



9e( x ) = g 



5 -13a; x(l - 9a;) , 
(1-x) 2 + (1-x) 3 



CoK,m 2 . m 3) 



m 



1 ln (™?//4) 



(mf — ml)(ml — m|) 



+ (mi <-> m 2 ) + (mi <-> m 3 ) 



/ 2 2 2\ ^ 1 

c 2 (m 1 ,m 2 ,m 3 ) = - - - 



m 



f ln(m2/^) 



mt 



m|)(mf — m| 



+ (mi <-> m 2 ) + (mi <-> m 3 ) 



do(m 1 ,m 2: m 3 , m 4 



m 2 ln(m 2 //i|) 



(mf — m|)(mf — m|)(mf — m 2 ) 



+ (mi <-> m 2 ) + (mi <-> m 3 ) + (mi <-> m 4 ) 



d 2 (ml,ml,ml,ml) 



m\ ln(m 2 //4) 



(ml — m|)(m 2 — m§)(m 2 — m 2 ) 



+ (mi <-> m 2 ) + (mi <-> m 3 ) + (m x <-> m 4 ) 



29 



APPENDIX D: FORM FACTORS 

In Tables | and [TT]we summarize the two different sets of form factors discussed 
which are related via 



F T {q 2 ) = -(M B + M P )s{q 2 ), 



V(q 2 ) = (M B + M v )g(q 2 ), 



f(q 2 ) 



Mb + M\ 



A 2 (q 2 ) = -(M B + M v )a + {q 2 ), 



A (q 2 ) 



, g 2 a_(g 2 ) + f(q 2 ) + (M| - M 2 )a + (q 2 ) 



2Mv 



Ti(q 2 ) = -\g + {q 2 ), 



T2{q 2 ) = ~- 



g+(q) + 



2, q 2 9-(q 2 ) 



M 2 B -Ml 



Uq 2 ) 



g-(q 2 ) 



2 , (M 2 -M 2 )h(q 2 ) 



30 



TABLE I. The B -» and B -> 7r(p) form factors of Melikhov and Nikitin g| 



Form factor 



B ^ K 

9 N "2-32 



7T 



-2.30 



s(q 2 ) 



0.36 1 



-0.30 1 



6.88 2 

„2 \ -2-27 



6.71 2 



0.06 GeV- 1 ( 1 



6.85 2 



-2.28 



0.29 1 



-0.26 1 



6.71 2 

„2 



6.55 2 



0.05 GeV 



-1 



6.68 2 



-2.31 



Form factor 



B -> K* 



2 \ -2.61 



B~ 



-2.75 



a + (q 2 ) 
a-(q 2 ) 
HQ 2 ) 
9+(q 2 ) 

g-(q 2 ) 



0.048 GeV" ] | 1 
1.61 GeV ^1 ■ 
-0.036 GeV" 



6.67 2 



q 2 q 4 S ~ V 
5.86 2 7.66 4 

2 \ -2.85 

7.332 J 

„% \ -2.72 



0.041 GeV" 1 ( 1 



0.0037 GeV" 2 ( 1 



-0.28 1 



6.98 2 

„2 \ "3.28 



6.572 

,2 \ -2'62 



0.24 1 



6.67 2 

2 \'-2.58 

T 
6.592 



0.036 GeV 



1.10 GeV 1 



-1 



6.552 



+ 



-1 



-0.026 GeV 



0.03 GeV" 



-1 



5.592 7 1Q 4 

2 \ -3.04 

1 



7.292 

„2 \ -2-85 



1 



6.88 2 

2 \ -3.42 

0.003 GeV -2 ( 1 - 



-0.20 1 



6.432 J 

,2 \ -2.76 



0.18 1 



6.57 2 

2 \'-2.73 

g 

6.50 2 
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TABLE II. The B -> JfW form factors of Colangelo ei a/. @, with M = 5 GeV. As for the 
B — > 7r(p) transition, we use the form factors listed below with M = 5.3 GeV in F\ and Ft- 



Form factor 



B -> K 



Fi(q 2 ) 
Fo(q 2 ) 
F T (q 2 ) 



0.25 1 



0.25 1 



-0.14 1 



M 2 



q 

M 2 

49 

-l 
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Form factor 



B -> K* 



V(q 2 ) 

Ax(q 2 ) 
A2(q 2 ) 

Ao(q 2 ) 

Ti(q 2 ) 

T2{q 2 ) 
Ts(q 2 ) 



0.47 1 



25 



0.37(1 - 0.023a 2 ) 
0.40(1 + 0.034a 2 ) 

0.30 1 



0.19 1 



4.8 2 

„2 



-1 



5.3 2 , 
0.19(1 - 0.02a 2 ) 

0.30(1 + o.oia 2 ) 
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